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04 | We report our first results of photometric and spectroscopic observations for Tycho's supernova rem- 

nant (SNR Tycho) to search for the companion star of a type la supernova (SN la). From photometric 
r~| | observations using Suprime-Cam on the Subaru Telescope, we have picked up stars brighter than 22 mag 

O |. (in V-band) for spectroscopy, which are located within a circular region with the radius of 30" around the 

center of SNR Tycho. If the ejecta of young supernova remnants, such as SNR Tycho, have a sufficient 
amount of Fe I, we should be able to detect absorption lines at 3720 A and 3860 A associated with transi- 
^ ' tions from the ground state of Fe I in the spectrum of the companion star. To identify the companion star 

of a SN la using these characteristic absorption lines of Fe I, we made optical low-resolution spectroscopy 
of their targets using FOCAS on the Subaru Telescope. In our spectroscopic observations, we obtained 
I \ spectra of 17 stars in the SNR Tycho region and compare them with template stellar spectra. We detect 

j> • significant absorption lines from two stars at 3720 A. In addition, spectra of four stars also have possible 

' absorption lines. Since widths of their absorption lines are broad, it is likely that the detected absorptions 

lO ■ are due to Fe I in the expanding ejecta of SNR Tycho. However, none of stars exhibits a clear red wing in 

' the observed profiles of the absorption, though a star in the background of the SNR should show it. Hence, 

we suggest another interpretation that the detected absorption lines might be caused by the peculiarity 
\»Q ■ of stars. Though Tycho(G) was raised as a candidate of the companion star of Tycho's supernova in a 

previous research, we do not obtain evidence supporting the previous research. On the contrary, a star 
named Tycho(E) has the absorption line at 3720 A and its projected position is close to the center of 
SNR Tycho. Based on our observations, Tycho(E) is a new candidate as the companion star of Tycho's 
supernova. 
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1. Introduction are the double degenerate model (DD) and the single 

degenerate model (SD). The DD scenario assumes that 

Type la supernovae (SNe la) have been used as stan- merging of double C+O white dwarfs with the combined 

dard distance indicators and have become a powerful tool mass exceeding the Chandrasekhar mass limit, which is 

for studying the expansion of the Universe, since they about 1.4 Mq, results in SN la explosion (e.g., Iben & 

exhibit similar maximum brightness and uniform light Tutukov 1984; Wcbbink 1984). However, the DD has 

curves. They are used to constrain the cosmological pa- been suggested to lead to an accretion-induced collapse 

rameters and provide the evidence for the existence of rather than a SN la (e.g., Nomoto & Iben 1985; Saio & 

dark energy (Schmidt et al. 1998;Perlmutter et al. 1999; Nomoto 1985). On the other hand, in the SD scenario, 

Knop et al. 2003; Astier et al. 2006; Riess et al. 2006). a white dwarf accretes the gas from the companion non- 

Despite of their importance in the cosmological stud- degenerate star in a binary system and increases its mass 

ies, the progenitor system of SNe la has not been fully up to the Chandrasekhar limit to explode as a SN (e.g., 

understood yet and is still debated (Branch et al. 1995). Nomoto et al. 1997; Han & Podsiadlowski 2004) . Hachisu, 

There are two models for the progenitor system, which Kato, &; Nomoto (1996) proposed a new plausible progen- 
itor system based on the optically thick wind theory of 



Based on data collected at Subaru Telescope, which is operated 
by the National Astronomical Observatory of Japan. 



mass-accreting white dwarfs. In this model, the accretion 
to the white dwarf at a rate exceeding the rate of hy- 
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drogen burning near the surface is inhibited by blowing a 
strong wind. In this way, the wind avoids the formation 
of a common envelope and the white dwarf can increase 
its mass up to the Chandrasckhar limit. They suggested 
two main evolutionary paths to SNe la, depending on the 
evolutionary stage of the companion star: a helium-rich 
supersoft X-ray source channel (Hachisu et al. 1999a) 
with the companion of a main sequence star and a wide 
symbiotic channel (Hachisu et al. 1999b) with a red giant. 

The above evolutionary scenarios for SNe la can be 
tested by observations. For example, a companion star 
should remain in the vicinity of the explosion site after 
the explosion of a SN la in the SD scenario. If the com- 
panion star is identified by observations, the SD scenario 
will be supported. 

Ruiz-Lapuente et al. (2004) (hereafter Ruiz04) reported 
results of their survey of the central region of SNR Tycho 
to search for the companion star. They measured veloc- 
ities and distances of stars in the vicinity of the center 
of the SNR and found a G2IV star moving at a velocity 
more than three times greater than the velocity dispersion 
of disk stars in that region. As a result, they claimed that 
this star was the companion star. (They called this star 
"Tycho(G)".) This type of companion stars contradicts 
the evolutionary paths proposed by Hachisu et al. (1999a, 
b) and seems to suggest an alternative evolutionary path 
to a SN la. Fuhrmann (2005) has argued that Tycho(G) 
might be a thick-disc star which is coincidentally passing 
the vicinity of the supernova remnant, though the prob- 
ability is extremely low. Therefore, the physical associa- 
tion of Tycho(G) with the SNR needs to be confirmed 
by other methods. To this end, we performed photo- 
metric and spectroscopic observations of stars including 
Tycho(G) with SUBARU telescope and present our re- 
sults in this paper. 

Ozaki & Shigeyama (2006) (hereafter OS06) proposed 
an alternative method to investigate the association of a 
star with the SN ejecta. When a star is in the vicinity of 
the center of the expanding ejecta, photons emitted from 
the star toward the observer are absorbed only by the 
ejecta moving toward the observer. Hence the absorption 
lines exhibit broad wings present only in the blue-shifted 
side, while a background star exhibits very broad wings 
in both of the blue and red sides of the absorption lines. 
The possibility that another irrelevant star happens to be 
inside the ejecta of SNR Tycho is only 0.1 % as estimated 
from the density of stars there. Therefore a star inside 
SNR Tycho is very likely to be the companion star of 
Tycho's SN. In SNR Tycho, atoms are expected to be in 
the ground states. Thus useful absorption lines in the 
optical wavelength range can be formed by Fe I at 3720 
A and 3860 A. 

Tycho's supernova (SN 1572) is one of the type la su- 
pcrnovae in our galaxy. The distance to SNR Tycho 
is rsj 3 kpc, derived from its X-ray spectrum (Fink et 
al. 1994). From the X-ray imaging observations with 
Chandra, the center of the remnant is determined to be (a, 
8) = (00 h 25 m 19?40, 64°08'13."98) (J2000.0) to minimize 
the ellipticity of the blast wave of the SNR (Warren et al. 



2005). They also estimated the diameter of SNR Tycho 
as 8'. SNR Tycho is one of a few well-studied young rem- 
nants with known age and the remnant is so young that 
it could retain enough amount of Fe I to form absorption 
lines in the spectrum of the companion star. 

In §2, we report two kinds of observations, imaging and 
spectroscopy. In §2.1, we describe results of imaging ob- 
servations and explain how we determine the region to 
search the companion star. In §2.2, we describe results 
of spectroscopic observations and data reduction. In §3, 
we explain the method to estimate the amount of absorp- 
tion lines of Fe I in terms of equivalent widths. In §4, 
we summarize the results of our observations. Finally, we 
discuss the detected absorption lines and the possibility 
of identification of the companion star of Tycho's SN in 
§5. 

2. Observations 

2.1. Photometric Observations 

2.1.1. Observations 

Stars in the region of SNR Tycho was observed with 
Suprime-Cam (Miyazaki et al. 2002) at the prime focus 
of 8.2 m Subaru telescope on Nov. 4, 2002. Suprime- 
Cam is equipped with ten MIT/LL 2048 x 4096 CCDs 
arranged in a 5 x 2 pattern to provide 34' x 27' FOV 
with a pixel size on the sky of 0.201" x 0.201". We used 
V and i? c filters of Johnson-Cousins systems and %' filter of 
SDSS (Fukugita et al. 1996). Three 10 sec and three 100 
sec exposures were taken using each filter. We divided 
exposure into three and dithered to remove cosmic rays 
and to avoid bad pixels. The individual dithered images 
were coadded. The seeing size was about 0.8" . 

2.1.2. Reduction and Photometry 

We analyze a 6' x 6' region around the center of SNR 
Tycho in images taken with Suprime-Cam. The size of 
the field is slightly larger than the field of view of FOCAS. 
The standard data reduction, including overscan subtrac- 
tion and flatfielding, is performed with reduction software 
nekosoft for Suprime-Cam (Yagi et al. 2002). We smooth 
images with Gaussian kernel and search peaks to iden- 
tify stellar positions. If the signal-to- noise ratio of a peak 
is higher than 5, then we identify the peak as an object. 
Photometry is then performed with a 2" </> aperture. Since 
we did not observe standard stars, the flux calibration of V 
and i? c -band is performed with the published photometry 
in SNR Tycho field (Ruiz04). For i'-band, we transform i' 
mag into I c mag of Johnson-Cousins systems by using the 
prescription in Fukugita et al. (1996). To calculate the 
zero point of I c mag, we use B and V mag of Ruiz04 and 
the color of B - V and V - I c of Soderblom et al. (1998). 
To avoid saturations and to increase signal-to-noise ratio 
as much as possible, the total exposure times are chosen 
to be 30 sec for bright objects with the magnitudes of 
16.0-18.5 mag (V), and 300 sec for faint objects with the 
magnitudes of 18.5-24.0 mag. The limiting magnitudes (5 
a) are 24.0 mag (V), 23.5 mag (R c ), and 23.5 mag (J c ) 
for 300 sec. We make a star catalog from these images. 
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Fig. 1. Color-magnitude diagram in V — I c versus V. m(V) 
is the apparent magnitude in V-band. Filled blue circles are 
stars in the candidate region. Filled green triangles are ob- 
served stars as candidates of background stars of SNR Tycho. 
We select blue stars as much as possible, because they are 
luminous. 




Fig. 2. Color-Color diagram in V — Rc versus V — I c - Filled 
blue circles are stars in the candidate region. Filled green 
triangles are observed stars as candidates of background stars 
of SNR Tycho. 

2.1.3. Results of Photometry 

The final catalog contains 904 stars with photome- 
try in V ,R Cl and I c bands. This catalog is limited by 
^-magnitudes. A color-apparent magnitude diagram of 
V — I c versus V is shown in Figure 1 and a color-color 
diagram of V — R c versus V — I c is shown in Figure 2. We 
show a y-band image taken with Suprime-Cam around 
the center of SNR Tycho in Figure 3 and a V^-band image 
which corresponds to the field of view of our spectroscopic 
observations in Figure 4. These observations are described 
in detail in §2.2. 

2.1. 4- Candidate region for the companion star 

We select candidates for the companion star of Tycho's 
SN using the projected distances from the center of SNR 
Tycho. We estimate the velocity of the companion star 
from the equation of motion of two bodies. In the SD 




Fig. 3. A V-band image around the center of SNR Tycho 
taken with Suprime-Cam for 300 sec exposure. The 30" ra- 
dius circle is the candidate region. The stars designated with 
capital alphabets are candidate stars of a companion star of 
Tycho's SN. The center of SNR Tycho is marked by a gray 
star symbol. 
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Fig. 4. A V-band image around the center of SNR Tycho 
with Suprime-Cam for 300 sec exposure. The area of the 
image is almost the same as the field of view of spectroscopic 
observations with FOCAS. The 1' diameter circle shows the 
candidate region for the companion star. The stars designated 
with lower case alphabets are spectroscopically observed to 
detect Fe I in the SNR. The center of SNR Tycho is marked 
by a star symbol. 
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model, there should have been a stellar binary system 
in which a white dwarf and the companion star revolve 
around their center of mass before SN la explosion. We 
calculate the orbital velocity of the companion star us- 
ing the final mass and the period presented in Hachisu & 
Kato (2001). We obtain the velocity of ~200 km/s or less. 
After the explosion, only the companion star remains and 
is supposed to move with the same velocity. Four hundred 
and thirty years have passed since the explosion of Tycho 
SN1572. The distance from the Earth to SNR Tycho is ~3 
kpc (§1). As a result, the companion star may have moved 
at most by ~10" from the center of SNR Tycho. Taking 
account of uncertainties in determination of the center of 
SNR Tycho, we regard a circular region with the radius 
of 30" as the candidate region. A y-band image of this 
region taken with Suprime-Cam is shown in Figure 3. We 
designate the stars in the candidate region with capital 
alphabets same as Ruiz04 except Tycho(Y), which has no 
data in Ruiz04 (see Figure 3). The stars out of the region 
are designated with lower case alphabets (see Figure 4). 
We select 8 stars brighter than 22 mag (V) in this region 
for spectroscopic observations. 

2.2. Spectroscopic Observations 

2.2.1. Observations 

We carried out optical spectroscopy of stars in the 
region of SNR Tycho using multi-object-spectroscopy 
(MOS) mode of the Faint Object Camera and 
Spectrograph (FOCAS) (Kashikawa et al. 2002) attached 
to Subaru Telescope. One mask of FOCAS/MOS covers a 
6' <j) aperture diameter field of view. The MOS slit lengths 
were 8."0 and slit widths were l."0. We used the 2nd order 
of the grism of 150 lines/mm and the L550 order-sorting 
filter. This setting gave wavelength coverage from 3600 A 
to 5200 A and the spectral resolution was R~400 at 3800 
A. The pixel scale of the CCDs is 0."1 /pixel. We binned 
3 pixels in the spatial direction. In order to get a higher 
throughput in short wavelength, observations were made 
without the atmospheric dispersion corrector (ADC). 

Observations were carried out on Sep. 14 and Sep. 15, 
2006. We observed 8 stars in the candidate regions within 
30" 4> radius and other 25 stars located outside of the can- 
didate region. The reason why we observed the outside 
stars is to confirm whether sufficient Fe I remains in SNR 
Tycho. The two brightest stars, (A) and (B) in the candi- 
date region, are not regarded as candidates because these 
stars arc so bright that they are very likely to be in the 
foreground of SNR Tycho. Seeing size during observations 
was ~0."6 - 0."8. The exposure time of each frame was 
0.5 hours and the total effective exposure time was 12.5 
hours. We dithered ±0."2 for each exposure to avoid bad 
pixels. We made two MOS masks to observe all of the 
candidate stars. The bright 5 targets including Tycho(G) 
were observed using Maskl. The exposure time of Maskl 
was 2.5 hours (0.5 hours x 5). Maskl was used at low 
elevation because targets of Maskl were relatively bright. 
The other 4 targets including Tycho(V) were observed us- 
ing Mask2. The exposure time of Mask2 was 10.0 hours 
(0.5 hours x 20). The other 24 targets were observed 



using both masks. 

2.2.2. Data and reduction 

We analyze obtained spectra using IRAF software 1 . 
Subtracting the bias level and flatficlding are carried out. 
Then, each spectrum is extracted from multislit images. 
Each spectrum of 0.5 hour exposure time is combined with 
imcombine of IRAF task. At the same time, this task re- 
jects cosmic rays and bad pixels above 3 a. We carry 
out skysubtraction to use a 2nd order of spline3 func- 
tion. We calibrate wavelengths using Th-Ar lamp. The 
one-dimensional spectra are made by apall of IRAF task 
and are flux-calibrated using spectra of a standard star, 
BD+28d4211, obtained with the long-slit spectroscopy 
mode of FOCAS with a 2."0 slit width. Then, we cor- 
rect these spectra for the Galactic extinction and slit loss. 
This correction is described in the next section in detail. 
We analyze 4 spectra Tycho(C), (E), (F) and (G) in the 
candidate region and spectra of 13 stars outside the can- 
didate region. The other spectra do not have good signal- 
to-noise ratios. 

Spectra of 21 stars before correction of the Galactic 
extinction and slit loss are available from the web page 
(http://www.ioa. s. u-tokyo.ac.jp/~iharayt/SNRTycho/). 
Out of 21 stars, 8 stars are in the candidate region and 
13 stars are outside the region. The above 21 spectra 
include 4 spectra of Tycho(D), (I), (V) and (Y) which 
do not have good signal-to-noise ratios (<3.0 at 3800 A) 
and we did not use them in this work. We analyzed the 
other 17 spectra with good signal-to-noise ratios (>3.0 
at 3800A). They are shown in Figure 6 after the flux 
correction. 

3. Method to detect Fe I absorption lines 

3.1. Calibration of observed spectra 

We need to take into account two causes of the flux 
loss at blue wavelengths. One of them is the reddening 
due to the Galactic extinction. It is difficult to accurately 
estimate the Galactic extinction because of unknown dis- 
tances to observed stars. The other is slit loss, which is 
expected to be a function of wavelength because we did 
not use ADC. We set the central wavelength of slit acquisi- 
tion at 3900-4100 A. The flux out of this wavelength range 
was gradually lost from the slit due to atmospheric dis- 
persions. Atmospheric dispersion becomes larger at short 
wavelengths compared with that at long wavelengths. We 
have to handle these two complex flux losses at the same 
time. Therefore, we make flux corrections by multiply- 
ing a following linear function so that a template stellar 
spectrum traces the observed spectrum. 

/(A) = l-ax(c-A) (1) 

where a is a variable and we correct the flux at shorter 
wavelengths than c. The values of a are listed up in Table. 

1 IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of 
Universities for Research in Astronomy, Inc., under cooperative 
agreement with the National Science Foundation 
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1. The value of c is fixed to 3840 A, because we avoid the 
wavelength range of 3840-3880 A where we expect Fe I 
absorption. 

Further details of spectral fitting and measuring possi- 
ble absorption lines are also described in §3.2 and §3.3. 

3.2. Spectral types of observed spectra 

We compare the observed spectra with the template 
spectra by Jacoby, Hunter, & Christian (1984) (hereafter 
Ja84). They incorporated spectra for many stars having 
spectral classes O-M and luminosity classes V, IE, and I 
into their library Their wavelength coverage is 3510-7400 
A and the resolution is ^4 A. Ja84 is the most suitable 
catalog compared with other spectral libraries, for exam- 
ple, Gunn & Stryker 1983, Silva & Cornell 1992, Soubiran, 
Katz & Cayrel 1998, or Pickles 1998, in terms of the num- 
ber of spectra, spectral resolution, and wavelength cover- 
ages. 

We determine spectral types by finding a template spec- 
trum that most closely matches the observed spectrum of 
each star in the wavelength range of 3900-4400 A. We have 
chosen this wavelength range because Fc I absorption lines 
due to SNR are expected at wavelengths shorter than 3900 
A, and also because the slit loss becomes large at around 
4400 A and increases at longer wavelengths. The res- 
olution of our observing mode was not good enough to 
distinguish whether the star was a dwarf or a giant. In 
other words, spectra of a dwarf and a giant are similar in 
low resolution. Thus we used only dwarf spectra for spec- 
tral fitting. This is sufficient for our aim to find possible 
absorption features of the SNR in observed spectra. 

We first compare the observed spectra with some tem- 
plates to see global line shapes at 3900-4400 A. Then we 
have chosen several template spectra by referring to the 
shapes of characteristic absorption lines such as the Ca 
H&K lines at 3968 A and 3933 A, the H<5 line at 4101 A, 
the CH line at 4300 A, and the H7 line at 4340 A. Finally, 
we calculate x 2 residuals of several candidate template 
spectra to determine the spectral types and the flux levels. 
To minimize possible errors in 3700-3900 A for detection 
of the Fc I absorptions, wc determine the spectral type 
by choosing the template spectrum with the minimum x 2 
residuals in the wavelength ranges of 3740-3840 A and 
3900-4100 A. In this x 2 fitting, the template spectra are 
binned to match the resolution to that of the observed 
spectra. The spectral types thus determined are listed in 
Table 1. In this fitting, we do not consider radial veloc- 
ities of stars because they are negligible compared with 
our spectral resolution. 

3.3. Measurement of absorption lines due to Fe I in the 
ejecta 

According to OS06, the spectrum of the companion star 
of Tycho's SN should show broad absorption lines of Fe I in 
the blue sides of 3720 A and 3860 A. In addition, the same 
Fe I in the ejecta should imprint broad absorption lines 
on both of the blue and red sides of 3720 A and 3860 A 
in the spectra of background stars of SNR Tycho. Figure 
5 shows the predicted profiles of these absorption lines 
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Fig. 5. The profiles of absorption lines in the spectrum of 
the companion star due to Fe I in the SN ejecta predicted for 
different values of the density of the ambient medium (Red 
solid line : n a = 2.0 cm -3 , Green dashed line: 1.8 cm -3 , Blue 
dotted line: 1.5 cm" 3 ). The companion star is assumed to be 
located within 10" from the center of SNR Tycho. 

of the companion star based on the W7 model (Nomoto, 
Thielemann, & Yokoi 1984) for several densities n a of the 
ambient medium. The profile is sensitive to n a of SNR 
Tycho. 

In this calculation, the hydrodynamical code is revised 
from OS06. We adopt the original distribution of elements 
in W7 instead of the simplified distribution used in OS06. 
As a result, elements such as Ca, Ar, S are added. Thus, 
we add new emission line data for the additional elements 
and revise line data of the rest elements. Moreover, the 
data set of collisional ionization cross section is updated. 
As we relax the assumption in OS06 that shock heated 
atoms are automatically ionized to some extent, more ion- 
izing photons are emitted from atoms in lower ionization 
stages. As a result, for n a =l. 5-2.0, which is higher than 
that calculated in OS06, the expected absorption line pro- 
files become shallower and weaker than those in OS06. 

We measure equivalent widths of possible absorption 
lines. This equivalent width is different from the usu- 
ally used equivalent width in the sense that the template 
spectrum is used for the reference instead of the contin- 
uum spectrum of the same star. Therefore, we introduce 
the term "Excess Absorption Measure" (EAM). The EAM 
here is defined as 

EAM =(U- /° bs(A) ... ) d\, (2) 

J \ template V A J / 

where the observed flux is denoted by -Fobs and the tem- 
plate flux ^template- We bin the spectra in order to reduce 
the photon noise and the readout noise. The binning size 
is 4 pixels (5.6 A per bin ). 

We consider two error sources. One is an error from 
the template fitting, and the other is a statistical error 
from photons and readout noise. The fitting error is es- 
timated by root mean square values (R.M.S.) (~Efa) of 
the residuals in the wavelength ranges where we expected 
no absorption lines, i.e., in 3650-3680 A, 3740-3840 A and 
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3880-3920 A. The photon noise is estimated by the num- 
ber of electrons originated from the stellar photons and 
the background photons at each wavelength, and signal- 
to-noise ratio (S/N) is given by the following formula, 



S/N 



-Est, a — 



VS + 37V + 30-2' 
1 



S/N' 



(3) 



(4) 



where S is the number of photons of an object per 3 pixels, 
N is the number of photons of the sky background per 
pixel averaged over adjacent 18 pixels, of which 9 pixels 
are on the right side of the signal and 9 pixels are on the 
left. er e is the readout noise, which is 4 e~ r.m.s/pixel for 
FOCAS. The gain is 2.1 e"/ADU. 

We combine these two error sources and the total error 
E a u is estimated by the following formula, 



-Bail - \l E Rt +E stli , 



(5) 



assuming that the residuals of the fitting follow the 
Gaussian distribution. We plot the ratios of E a \\ to F b s 
per bin in Figure 6-22. 

We calculate EAMs in five wavelength ranges, which are 
the blue and red sides of 3720 A and 3860 A , and around 
3800 A. The EAMs are calculated by integrating 4 bins 
(22 A) which correspond to the absorption widths pre- 
dicted for the blue wing of the companion star by OS06. 
The errors are also accumulated in the same bins. We 
summarize the results with errors in Table. 1. 

4. Results 

4-1. Spectral types 

The results are summarized in Figures 6-22 and Table 
1. In Table 1, we list names of stars, fitted spectral types, 
the number of masks used in the spectroscopy, signal-to- 
noise ratios, fitting parameters (% 2 and a), magnitudes in 
V-band, projected distances from the center of the SNR, 
and EAMs (of the blue and red wings of lines at 3720 
A, at 3860 A, and 3790-3810 A). The EAM of 3790-3810 
A (labeled as 3800 in Table 1) is used to see how the 
template spectra match the observed spectra. Although 
we take account of only main sequence stars in the spec- 
tral fitting, spectral types from this method are differ- 
ent from the results of Ruiz04. For example, we classify 
Tycho(G) as a F8V star while Ruiz04 classified it as a 
G2IV star. Furthermore, we classify Tycho(E) as a F8V 
star, while Ruiz04 classified it as a K2-3ITI star, and we 
classify Tycho(F) as a F4V star, while Ruiz04 classified 
it as a F9ITI star. Ruiz04 classified stellar types using the 
spectra at around 6500 A. On the other hand, we obtained 
spectra with lower resolutions in the range of 3900-4400 
A (see §3.2). 

4-2. Absorption feature 

4.2.1. Tycho(G) 

Since Tycho(G) was raised as the companion star of 
Tycho's SN by Ruiz04, we discuss this star first. As 



shown in Figures 9 and Table 1, the observed spectrum 
of Tycho(G) does not exhibit any significant broad ab- 
sorption lines (EAM(3720B)=0.59±0.53A) expected due 
to Fe I in the ejecta and is consistent with the template 
spectrum of a F8V star. Therefore we do not obtain sup- 
porting evidence for Tycho(G) as the companion star of 
Tycho's SN. Nevertheless, this result does not refute the 
association of this star with Tycho SNR. There might re- 
main only a fraction of Fe I in the ejecta insufficient to 
make broad absorption lines in the spectrum of the com- 
panion star. 
4.2.2. Tycho(E) 

Significant broad absorption lines with positive EAMs 
are detected from stars Tycho (E), (aa), (ad), (ah), (ai), 
and (ak). Since only Tycho(E) is a candidate star for the 
companion of Tycho's SN in terms of its position, we take 
this star as an example and discuss the absorption profile 
of the spectrum. First, we discuss the spectrum around 
3720 A. There is a significant absorption with the EAM 
of 3.0 A (5.5 cr) in the blue side of 3720 A. The profile of 
the absorption resembles the theoretical prediction with 
n =1.8 enr 3 in Figure 5. The absorption line appears 
to extend to 3685 A. On the other hand, no significant 
absorption is detected around 3860 A. If the absorption 
at 3720 A is due to Fe I, the absorption with about a 
half of the EAM at 3720 A is expected at 3860 A. Since 
the accuracy of fitting around 3860 A is rather poor, we 
also fit the observed spectrum with spectra of F7V and 
F9V stars in addition to F8V (Figure 23) to examine the 
influence of the adopted spectral type upon detections of 
broad absorption lines. As shown in Figure 23, the shapes 
of the three spectra at around 3860 A differ from each 
other, while the shapes at around 3720 A are similar. If 
Tycho(E) had a spectrum of a F7V star as its original 
spectrum, the EAM in the blue side of 3860 A would be- 
come large (~2.5 A). However, the observed spectrum is 
fitted with F8V better than F7V. If we fitted with the 
middle type between F7V and F8V, we may detect some 
absorption at 3860 A. Because of a large variation of stel- 
lar spectra at around 3860 A, it is difficult to quantita- 
tively discuss the relatively weak absorption line at 3860 

A. 

Ruiz04 classified this star as K2-3IH and estimated the 
distance greater than 20 kpc. This apparently contra- 
dicts our classification based on the spectrum observed in 
3900-4400 A. The spectrum of a catalogued star of K2II 

I (or K3HI) is quite different from that of F8V not only 
in this wavelength range but also at around 6500 A that 
Ruiz04 used to deduce the spectral type. As shown in 
the color-color diagram of the observed stars (Fig. 2), 
the colors of Tycho(E) are V -R c = 0.9 and V-I c = 1.9, 
similar to that of Tycho(G). Both of the stars are classi- 
fied as F8V. The colors of Tycho(C), which is classified 
as M1V, are V - R c = 1.2 and V - I c = 2.4. If a K2-3I 

II star is located at the distance of 20kpc, the colors of 
K2-3HI should be V - R c > 1.2 and V - I c > 2.5 because 
this star is affected by larger the Galactic extinction than 
other stars which are located at the distance of 3kpc or 
less. As shown in the color-color diagram, stars are lo- 



Table 1. Summary of observations 



Star name 


Type 


V mag 


Distance a 


Mask b 


S/N c 


x 2d 


a d 


EAM(3720B) 


EAM(3720R) 


EAM(3800) 


EAM(3860B) 


EAM(3860R) 


TvchofG) 


M1V 


18.9 


4.5" 


1&2 


7.5 


6.33 


0.0004 


52±0 96 


-1 08±0 92 

_L • KJ ^7 — 1 — \J • \J t-i 


37±0 83 

KJ ■ KJ 1 — 1 — KJ • KJKJ 


37±0 78 

KJ ■ KJ 1 — 1 — v7 ■ 1 v7 


75±0 77 

V7 • J KJ — 1 — Kt ■ 1 1 


TvchofD) 




20.6 


5.6" 


1&2 


2.8 
















TvchofE) 


F8V 


19.7 


10.8" 


1&2 


13.3 


7.45 


0006 

kj •kjkjkj kj 


3 00±0 55 

KJ ■ \J \J — 1 — V 7 • • ./ '.J 


56±0 52 

V / . kJ \J — 1 — KJ ■ \J i-t 


0.54±0.44 


-0.10±0.42 


-0 30±0 41 

V 7 ■ KJ KJ — 1 — V7 ■ I 1 


Tycho(F) 


F4V 


17.7 


26.2" 


1 


27.1 


26.0 


0.0023 


0.78±0.42 


0.09±0.41 


0.47±0.40 


0.13±0.39 


-0.27±0.39 


TvchoCG) 


F8V 


18.7 


26.2" 


1 


14.9 


11.9 


0020 


59±0 53 


-0 25±0 51 

V 7 . <—t -.7 — 1 — KJ • KJ J- 


-0 36±0 45 

KJ *KJKJ — 1 — KJ • 1 KJ 


-0 06±0 43 

KJ • KJ KJ — J — V ' • 1 f / 


-0 58±0 42 

KJ ■ KJ KJ — 1 — V ' > I 


Tvcho(T) 




21.8 


22.7" 


1&2 


2.3 
















Tycho(V) 




21.4 


24.8" 


2 


2.4 
















Tycho(Y) 




22.1 


24.2" 


1&2 


2.0 
















Tycho(aa) 


F8V 


20.3 


1.20' 


2 


6.7 


3.26 


0.0028 


3.06±1.07 


0.26±0.96 


0.51±0.71 


0.14±0.61 


-0.61±0.58 


Tycho(ab) 


G6V 


20.8 


1.26' 


1&2 


3.1 


1.39 


0.0007 


0.26±1.71 


1.24±1.58 


1.18±1.25 


0.05±1.05 


1.15±0.97 


Tycho(ac) 


F5V 


19.5 


2.40' 


1 


10.1 


6.97 


0.0023 


0.61±0.86 


-1.11±0.80 


-0.04±0.67 


0.82±0.62 


0.53±0.60 


Tycho(ad) 


G3V 


20.9 


1.27' 


2 


3.9 


1.45 


0.0004 


2.29±1.28 


1.69±1.20 


0.06±1.03 


-0.24±0.93 


0.24±0.89 


Tycho(ae) 


G6V 


20.1 


2.01' 


1&2 


6.7 


2.06 


0.0012 


-0.24±0.86 


0.67±0.82 


0.21±0.72 


0.94±0.67 


0.89±0.65 


Tycho(af) 


G2V 


20.3 


2.20' 


1&2 


7.3 


3.36 


0.0004 


-0.91±1.11 


-0.53±1.01 


0.80±0.83 


-0.28±0.77 


0.18±0.75 


Tycho(ag) 


F6V 


19.6 


1.00' 


1&2 


16.3 


6.47 


0.0015 


0.20±0.41 


-0.74±0.39 


0.17±0.34 


0.56±0.31 


0.02±0.30 


Tycho(ah) 


F8V 


19.7 


1.94' 


1&2 


13.2 


8.92 


0.0022 


2.91±0.42 


0.72±0.39 


0.10±0.33 


-0.01±0.29 


-0.53±0.28 


Tycho(ai) 


G3V 


21.0 


1.42' 


1&2 


3.8 


2.06 


0.0016 


3.86±1.94 


2.22±1.81 


-0.45±1.56 


2.06±1.49 


1.86±1.47 


Tycho(aj) 


F8V 


20.3 


2.38' 


1&2 


8.8 


4.53 


0.0023 


-0.58±0.84 


0.01±0.79 


-0.06±0.70 


0.31±0.66 


-0.57±0.65 


Tycho(ak) 


G3V 


20.9 


1.85' 


1&2 


3.9 


1.16 


0.0004 


2.67±1.60 


1.85±1.50 


-1.21±1.24 


1.14±1.11 


1.76±1.06 


Tycho(al) 


G6V 


20.1 


2.27' 


1&2 


6.8 


2.26 


0.0008 


0.52±0.97 


0.62±0.91 


-0.21±0.77 


0.13±0.70 


-0.13±0.68 


Tycho(am) 


F7V 


19.9 


3.74' 


1&2 


12.3 


6.56 


0.0007 


0.54±0.64 


-0.03±0.60 


-0.37±0.54 


-0.10±0.50 


0.62±0.49 



a The distance is the projected distance from the center of SNR Tycho which is {a, S) = (00 h 25 m 19?40, 64°08'13"98) (J2000.0). 
b We made two MOS masks according to stellar positions (see §2.2.1). S/Ns are measured at 3800 A per 1 pixel (1.4 A). 

d They are fitting parameters. Their details are written in §3.1 and §3.2. 

* The stellar types of Tycho(D), (I), (V) and (Y) are unknown for not good signal-to- noise ratios. 
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Fig. 6a. Tycho(C) Fig. 6b. Tycho(C), ratio 

Fig. 6. Comparisons of the observed spectra with template stellar spectra (Jacoby, Hunter, & Christian 1984). The left panels 
are our observed spectra (Red solid lines) and the fitted template spectra (Green dashed lines). The right panels show the ratio of 
observed spectra to the fitted template spectra. The dashed lines show lcr error per bin estimated by (5). According to Ozaki & 
Shigeyama (2006), absorption lines of Fe I are expected at ~3720 A and ~3860 A if stars are within or behind SNR Tycho, and 
those wavelengths arc shown by two dotted vertical lines. 
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Fig. 7. Same as Figure 6, but for Tycho(E). 
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Fig. 8b. Tycho(F), ratio 



Fig. 8. Same as Figure 6, but for Tycho(F). 
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Fig. 9a. Tycho(G) 
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Fig. 9. Same as Figure 6, but for Tycho(G). 





Fig. 11. Same as Figure 6, but for Tycho(ab). 
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Fig. 12a. Tycho(ac) 
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Fig. 12. Same as Figure 6, but for Tycho(ac). 
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Fig. 13a. Tycho(ad) 
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Fig. 13b. Tycho(ad), ratio 



Fig. 13. Same as Figure 6, but for Tycho(ad). 




Fig. 14. Same as Figure 6, but for Tycho(ae). 
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Fig. 15a. Tycho(af) 
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Fig. 15. Same as Figure 6, but for Tycho(af). 





Fig. 17. Same as Figure 6, but for Tycho(ah). 
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Fig. 18a. Tycho(ai) 
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Fig. 18. Same as Figure 6, but for Tycho(ai). 





Fig. 20. Same as Figure 6, but for Tycho(ak). 



No. 



Searching for a companion of SNIa 



13 



1.4 
1.2 
1 

0.8 

0.6 
0.4 
0.2 




3720 



3860 



3720 



3860 



- 






Tycho(al) — 
G6V 


1.2 
1 

0.8 


Ml 








0.6 










0.4 










0.2 



3700 3800 



3900 4000 
Wavelength (A) 



4100 



4200 



Fig. 21a. Tycho(al) 
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Fig. 21b. Tycho(al), ratio 



Fig. 21. Same as Figure 6, but for Tycho(al). 
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4100 



4200 



1.2 

1 

0.8 
0.6 
0.4 
0.2 





3650 





Tycho(am 




- 




F7\, 









































3700 



3750 3800 
Wavelength (A) 

Fig. 22b. Tycho(am), ratio 
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Fig. 22. Same as Figure 6, but for Tycho(am) 



catcd on a sequence that has some scatter in the region 
with V - R c > 3. Tycho(E) and Tycho(G) are located left 
on the sequence while the K2-3III star should be located 
right on the sequence because a M1V star Tycho(C) is lo- 
cated just on the sequence. Therefore, the classification of 
Tycho(E) by Ruiz04 contradicts the colors of Tycho(E). It 
is not at all impossible that Tycho(E) might have a quite 
peculiar spectrum if taking these results at face value. 
Further spectroscopic observations covering whole of the 
above wavelengths are needed to decide the spectral type 
of this star. 
4.2.3. Tycho(ah) 

The absorption lines are seen at 3720 A from the spec- 
tra of other stars, but most of their profiles arc different 
from that of Tycho(E). Tycho(ah) has a broad absorption 
line with the profile very similar to that of Tycho(E). In 
fact, the overall spectra of these two stars are quite similar. 
The projected distance of Tycho(ah) from the geometrical 
center of SNR Tycho is -2'. Thus if Tycho(E) is close to 
the SN site, then Tycho(ah) should have no physical in- 
teraction with the SN and vice versa. If Tycho(ah) is one 
of the background stars of the SNR, the spectrum should 



have a broad absorption line in the red side of ~3720 A in 
addition to the blue side. Our observations indicate that 
this would not be the case. Hence, the observed spectra 
of Tycho(E) and (ah) may indicate that they belong to 
a peculiar stellar type of F8V that has some deficit at 
around -3720 A compared with the F8V of Ja84. This 
deficit might be caused by metal richness of the star as 
the spectra of F8V stars with different metallicities in the 
catalogue of Pickels (1998) shown in Figure 24. We also 
compare with the F8V stars in the catalogue of Bagnulo 
et al (2003) shown in Figure 25. A lot of metal absorp- 
tion lines from 3720 A to 3900 A changing the shape of 
the spectrum in this region make it difficult to deduce the 
accurate EAMs of the blue wing of 3720 A and 3860 A 
absorption. To do this, we need to observe various nearby 
stars with the spectral type of F8V toward the direction 
of SNR Tycho and out of the region of SNR Tycho. 

5. Discussion 

We present a statistical study of the broad absorption 
lines at around 3720 A . The EAM distribution of the ob- 
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Fig. 24. Comparison of the template spectra of F8V stars 
with different metallicities from Pickels 1998. The red line 
shows the spectrum of a metal-weak F8V star, the green line 
shows a F8V star with the solar metallicity, and the blue line 
shows a metal-rich F8V star. The green line almost follows 
the spectrum of the F8V star of Ja84. The difference between 
the metal-rich and the normal star does not affect the amount 
of the absorption at ~3720 A. If we use the metal-weak star, 
the amount of the absorption is overestimated. 



Fig. 23. Comparisons of the observed spectrum of Tycho(E) 
(red solid curve) with the template spectra of F7V (green 
dashed curve), F8V (blue dashed curve), and F9V (purple 
dashed curve) stars of Ja84. The upper panel shows spectra 
in the whole observed range of wavelength. The bottom panel 
is a close-up of the spectra in a limited wavelength range 
3650-3900 A. The three template stars show similar behavior 
at around 3720 A, while they are significantly different at 
around 3860 A. Besides, a large difference is not seen at the 
other wavelengths. 

served stars is plotted in Figure 26 and the EAMs in the 
wavelength ranges of 3700-3720 A and 3790-3810 A are 
plotted with their magnitudes in V-band in Figure 27. 

In 3700-3720 A a broad absorption line due to Fe I 
is expected not only in the spectrum of the companion 
star but also in the spectra of the background stars, while 
in 3790-3810 A no absorption is expected. In 3790-3810 
A, the obtained EAM distribution is consistent with a 
scatter due to measurement errors around EAM -0 A 
(Figure 26b and 27b). On the other hand, there seem to 
be several stars which show positive EAMs of around 3 A 
in addition to stars without showing finite EAMs (Figure 
26a and 27a). 

Next, we discuss the profiles of absorption lines. The 
profile in the spectrum of Tycho(E) was discussed in §4. 
If Tycho(E) is the companion star, then the other stars 
with broad absorption lines should be located in the back- 
ground of SNR Tycho. These stars are expected to show 
broad absorption lines due to Fe I in the ejecta both in 
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Fig. 25. Comparison of the template spectra of F8V stars in 
Bagnulo ct al (2003). These spectra were obtained by mem- 
bers of the ESO Paranal Science Operations Team in order to 
make a library of high-resolution spectra of stars across the 
Hcrtzsprung-Russell diagram. The spectra by Bagnulo et al 
(2003) are converted into lower resolution spectra suitable for 
Ja84. We plot the four F8V stars in the catalog and one F8V 
star in Ja84 (red line). None of their spectra can exhibit the 
absorption feature like Tycho(E) and (ah). 
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-2 2 

Excess absorption measure (A) 

Fig. 26a. Distribution of EAM at 3700-3720 A 
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Fig. 26b. Distribution of EAM at 3790-3810 



Fig. 26. The EAMs versus the number of stars. There are two peaks at ~0.5A and ~3.0A at 3700-3720 A. On the other hand, 
there is only one peak at OA for 3790-3810 A. 




-4-2 2 

Excess absorption measure (A) 

Fig. 27a. Distribution of EAM at 3700-3720 
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Fig. 27b. Distribution of EAM at 3790-3810 A 



Fig. 27. The EAMs versus the magnitude in V-band. The dashed vertical lines indicate the mean value given by Figure 26. The 
stars of about 21 mag have large errors. 




-4 -2 2 4 

Excess absorption measure (A) 

Fig. 28a. Distribution of EAM at 3720-3740 A 
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Fig. 28b. EAM versus V-mag 



Fig. 28. Distribution of EAMs at 3720-3740 A same as Figure 26 and 27. We expect to detect the absorption lines from the 
background stars, but we cannot detect them clearly like the EAMs at 3700-3720 A. 
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the blue and red sides of the original wavelength 3720 
A. However, they exhibit very weak absorptions in the 
red side compared with that in the blue. In Figure 28, we 
show the EAM distribution in the range 3720-3740 A sim- 
ilar to Figures 26 and 27. There seem to be 4 stars having 
weak absorptions, but they are faint and the EAMs have 
large errors. In contrast, the other stars are distributed 
around A. 

There are 5 stars outside of the candidate region from 
which broad absorption lines are detected. Tycho(ah) has 
the best signal-to-noise ratio and is already discussed in 
§4.2.3. Tycho(aa) exhibits a similar spectrum but with 
smaller signal-to-noise ratios. Tycho(ad), (ai), and (ak) 
seem to have broad absorption lines with larger EAMs in 
the red side of 3720 A than Tycho(ah). These EAMs of 
Tycho(ad), (ai), and (ak) are about a half of the EAMs 
of the corresponding blue side absorption lines. However, 
large errors involved in these EAMs present us from de- 
ducing a clear conclusion. 

If positive EAMs are not caused by spectral features 
of peculiar stars, absorptions of Fe I is the most proba- 
ble interpretation. However, the observed absorptions are 
not exactly the same as predicted by the model of OS06. 
No red wings are significantly detected, while a significant 
blue wing is found for Tycho(ah) although its projected 
location is far from the center of SNR. Hence, if the origin 
of this absorption is Fe I, refining OS06 would be neces- 
sary. Their model assumes spherical symmetry distribu- 
tion of the ejected iron. But if distribution of neutral iron 
is clumpy, the Fe I absorption might not be able to be 
detected according to a location of a star because of lack 
of the iron in the line of sight. Also absorption profiles 
would be complicated. 

Finally, we discuss Tycho(E) as candidates of the com- 
panion star of Tycho's SN. Although there is no clear ex- 
ample which shows significant broad absorption lines both 
in blue and red sides of 3720 A which are expected for the 
background of the SNR, we confirm the blue-shifted broad 
absorption line at 3720 A only for Tycho(E) in the can- 
didate region. An absorption line with this unique shape 
can be formed by Fe I in the SN ejecta when the star is 
inside SNR Tycho. If the supernova explosion occurred at 
the center of SNR Tycho defined by Warren et al. (2005), 
the projected distance of Tycho(E) is considerably short 
from the center of SNR Tycho (-10"), while Tycho(G) 
is located away from the center of SNR Tycho (~25"). 
Ruiz04 suggested that the SN site was about 25" off the 
geometrical center of the SNR. It is not clear whether this 
SN site is consistent with the observed shape of the SNR in 
the context of hydrodynamical evolution of the blast wave 
in the ambient inhomogeneous medium. Therefore, it is 
likely that Tycho(E) is another candidate for the compan- 
ion star, considering the position from the center of SNR 
Tycho. 

6. Summary 

We have performed photometric and spectroscopic ob- 
servations to search for the companion star remaining in 



the remnant of a SN la (SNR Tycho). First, we selected 
candidate stars for the spectroscopic observations from 
photometric observations for the SNR Tycho region. Then 
wc obtained optical spectra of 4 stars in the candidate re- 
gion and 13 other stars outside of the candidate region. 
To identify the companion star of Tycho's SN, we com- 
pare the observed spectra in SNR Tycho with spectra of 
Ja84 at 3720 A, 3800 A, and 3860 A. If the position of a 
star is inside the ejecta of SNe la or in the background of 
the ejecta. absorption lines due to Fe I in the expanding 
ejecta may be seen. 

We are able to detect significant broad absorption lines 
at 3720 A for Tycho(E) and (ah). As the distribution 
of EAMs at 3720 A shows, spectra of Tycho(aa), (ad), 
(ai) and (ak) have possible absorption lines. Tycho(E) 
is located at the candidate region, while Tycho(aa), (ad) 
,(ah), (ai) and (ak) are located out of the region. However, 
we cannot detect the absorption at 3860 A at a sufficient 
significance level. It is difficult to detect them at 3860 A 
because of uncertainty of template stellar spectra at 3860 

A. 

The detected absorption lines in the blue side at 3720 
A are consistent with Fe I absorption of SNR Tycho. On 
the other hand, we do not obtain a clear broad absorption 
line in the red side of the original wavelength from a star 
in the background of the SNR. We cannot conclude if the 
detected broad absorption lines are due to Fe I in the SN 
ejecta. Another interpretation is that their absorption 
lines might be caused by the peculiarity of stars. 

We do not detect any broad absorption lines in the spec- 
trum of Tycho(G). Thus we do not obtain supporting ev- 
idence for the result of Ruiz04 from our observations. In 
the candidate region, Tycho(E) is the only star within the 
candidate region that shows a significant broad absorp- 
tion line. The spectral type of Tycho(E) determined by 
our observations and by Ruiz04 are apparently inconsis- 
tent. Our observations indicate that Tycho(E) is a new 
candidate as a companion star in the sense that this star 
is located close to the projected center of the SNR and 
shows a broad absorption line. 
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